The Concise Guide to PHARMACOLOGY 2023/24: Ion channels.

The Concise Guide to PHARMACOLOGY 2023/24 is the sixth in this series of biennial publications. The Concise Guide provides concise overviews, mostly in tabular format, of the key properties of approximately 1800 drug targets, and over 6000 interactions with about 3900 ligands. There is an emphasis on selective pharmacology (where available), plus links to the open access knowledgebase source of drug targets and their ligands (https://www.guidetopharmacology.org/), which provides more detailed views of target and ligand properties. Although the Concise Guide constitutes almost 500 pages, the material presented is substantially reduced compared to information and links presented on the website. It provides a permanent, citable, point-in-time record that will survive database updates. The full contents of this section can be found at http://onlinelibrary.wiley.com/doi/10.1111/bph.16178. Ion channels are one of the six major pharmacological targets into which the Guide is divided, with the others being: G protein-coupled receptors, nuclear hormone receptors, catalytic receptors, enzymes and transporters. These are presented with nomenclature guidance and summary information on the best available pharmacological tools, alongside key references and suggestions for further reading. The landscape format of the Concise Guide is designed to facilitate comparison of related targets from material contemporary to mid-2023, and supersedes data presented in the 2021/22, 2019/20, 2017/18, 2015/16 and 2013/14 Concise Guides and previous Guides to Receptors and Channels. It is produced in close conjunction with the Nomenclature and Standards Committee of the International Union of Basic and Clinical Pharmacology (NC-IUPHAR), therefore, providing official IUPHAR classification and nomenclature for human drug targets, where appropriate.

Implications of a temperature-dependent heat capacity for temperature-gated ion channels.

Temperature influences dynamics and state-equilibrium distributions in all molecular processes, and only a relatively narrow range of temperatures is compatible with life-organisms must avoid temperature extremes that can cause physical damage or metabolic disruption. Animals evolved a set of sensory ion channels, many of them in the family of transient receptor potential cation channels that detect biologically relevant changes in temperature with remarkable sensitivity. Depending on the specific ion channel, heating or cooling elicits conformational changes in the channel to enable the flow of cations into sensory neurons, giving rise to electrical signaling and sensory perception. The molecular mechanisms responsible for the heightened temperature-sensitivity in these ion channels, as well as the molecular adaptations that make each channel specifically heat- or cold-activated, are largely unknown. It has been hypothesized that a heat capacity difference (ΔCp) between two conformational states of these biological thermosensors can drive their temperature-sensitivity, but no experimental measurements of ΔCp have been achieved for these channel proteins. Contrary to the general assumption that the ΔCp is constant, measurements from soluble proteins indicate that the ΔCp is likely to be a function of temperature. By investigating the theoretical consequences for a linearly temperature-dependent ΔCp on the open-closed equilibrium of an ion channel, we uncover a range of possible channel behaviors that are consistent with experimental measurements of channel activity and that extend beyond what had been generally assumed to be possible for a simple two-state model, challenging long-held assumptions about ion channel gating models at equilibrium.

Calcium dependence of both lobes of calmodulin is involved in binding to a cytoplasmic domain of SK channels.

KCa2.1-3 Ca2+-activated K+-channels (SK) require calmodulin to gate in response to cellular Ca2+. A model for SK gating proposes that the N-terminal domain (N-lobe) of calmodulin is required for activation, but an immobile C-terminal domain (C-lobe) has constitutive, Ca2+-independent binding. Although structures support a domain-driven hypothesis of SK gate activation by calmodulin, only a partial understanding is possible without measuring both channel activity and protein binding. We measured SK2 (KCa2.2) activity using inside-out patch recordings. Currents from calmodulin-disrupted SK2 channels can be restored with exogenously applied calmodulin. We find that SK2 activity only approaches full activation with full-length calmodulin with both an N- and a C-lobe. We measured calmodulin binding to a C-terminal SK peptide (SKp) using both composition-gradient multi-angle light-scattering and tryptophan emission spectra. Isolated lobes bind to SKp with high affinity, but isolated lobes do not rescue SK2 activity. Consistent with earlier models, N-lobe binding to SKp is stronger in Ca2+, and C-lobe-binding affinity is strong independent of Ca2+. However, a native tryptophan in SKp is sensitive to Ca2+ binding to both the N- and C-lobes of calmodulin at Ca2+ concentrations that activate SK2, demonstrating that the C-lobe interaction with SKp changes with Ca2+. Our peptide-binding data and electrophysiology show that SK gating models need deeper scrutiny. We suggest that the Ca2+-dependent associations of both lobes of calmodulin to SKp are crucial events during gating. Additional investigations are necessary to complete a mechanistic gating model consistent with binding, physiology, and structure.

THE CONCISE GUIDE TO PHARMACOLOGY 2021/22: Ion channels.

The Concise Guide to PHARMACOLOGY 2021/22 is the fifth in this series of biennial publications. The Concise Guide provides concise overviews, mostly in tabular format, of the key properties of nearly 1900 human drug targets with an emphasis on selective pharmacology (where available), plus links to the open access knowledgebase source of drug targets and their ligands (www.guidetopharmacology.org), which provides more detailed views of target and ligand properties. Although the Concise Guide constitutes over 500 pages, the material presented is substantially reduced compared to information and links presented on the website. It provides a permanent, citable, point-in-time record that will survive database updates. The full contents of this section can be found at http://onlinelibrary.wiley.com/doi/bph.15539. Ion channels are one of the six major pharmacological targets into which the Guide is divided, with the others being: G protein-coupled receptors, nuclear hormone receptors, catalytic receptors, enzymes and transporters. These are presented with nomenclature guidance and summary information on the best available pharmacological tools, alongside key references and suggestions for further reading. The landscape format of the Concise Guide is designed to facilitate comparison of related targets from material contemporary to mid-2021, and supersedes data presented in the 2019/20, 2017/18, 2015/16 and 2013/14 Concise Guides and previous Guides to Receptors and Channels. It is produced in close conjunction with the Nomenclature and Standards Committee of the International Union of Basic and Clinical Pharmacology (NC-IUPHAR), therefore, providing official IUPHAR classification and nomenclature for human drug targets, where appropriate.

Structural basis of cytoplasmic NaV1.5 and NaV1.4 regulation.

Voltage-gated sodium channels (NaVs) are membrane proteins responsible for the rapid upstroke of the action potential in excitable cells. There are nine human voltage-sensitive NaV1 isoforms that, in addition to their sequence differences, differ in tissue distribution and specific function. This review focuses on isoforms NaV1.4 and NaV1.5, which are primarily expressed in skeletal and cardiac muscle cells, respectively. The determination of the structures of several eukaryotic NaVs by single-particle cryo-electron microscopy (cryo-EM) has brought new perspective to the study of the channels. Alignment of the cryo-EM structure of the transmembrane channel pore with x-ray crystallographic structures of the cytoplasmic domains illustrates the complementary nature of the techniques and highlights the intricate cellular mechanisms that modulate these channels. Here, we review structural insights into the cytoplasmic C-terminal regulation of NaV1.4 and NaV1.5 with special attention to Ca2+ sensing by calmodulin, implications for disease, and putative channel dimerization.

Ancestral reconstruction of protein interaction networks.

The molecular and cellular basis of novelty is an active area of research in evolutionary biology. Until very recently, the vast majority of cellular phenomena were so difficult to sample that cross-species studies of biochemistry were rare and comparative analysis at the level of biochemical systems was almost impossible. Recent advances in systems biology are changing what is possible, however, and comparative phylogenetic methods that can handle this new data are wanted. Here, we introduce the term "phylogenetic latent variable models" (PLVMs, pronounced "plums") for a class of models that has recently been used to infer the evolution of cellular states from systems-level molecular data, and develop a new parameterization and fitting strategy that is useful for comparative inference of biochemical networks. We deploy this new framework to infer the ancestral states and evolutionary dynamics of protein-interaction networks by analyzing >16,000 predominantly metazoan co-fractionation and affinity-purification mass spectrometry experiments. Based on these data, we estimate ancestral interactions across unikonts, broadly recovering protein complexes involved in translation, transcription, proteostasis, transport, and membrane trafficking. Using these results, we predict an ancient core of the Commander complex made up of CCDC22, CCDC93, C16orf62, and DSCR3, with more recent additions of COMMD-containing proteins in tetrapods. We also use simulations to develop model fitting strategies and discuss future model developments.

Non-Additive Effects of Binding Site Mutations in Calmodulin.

Despite decades of research on ion-sensing proteins, gaps persist in the understanding of ion binding affinity and selectivity even in well-studied proteins such as calmodulin. Site-directed mutagenesis is a powerful and popular tool for addressing outstanding questions about biological ion binding and is employed to selectively deactivate binding sites and insert chromophores at advantageous positions within ion binding structures. However, even apparently nonperturbative mutations can distort the binding dynamics they are employed to measure. We use Fourier transform infrared (FTIR) and ultrafast two-dimensional infrared (2D IR) spectroscopy of the carboxylate asymmetric stretching mode in calmodulin as a mutation- and label-independent probe of the conformational perturbations induced in calmodulin's binding sites by two classes of mutation, tryptophan insertion and carboxylate side-chain deletion, commonly used to study ion binding in proteins. Our results show that these mutations not only affect ion binding but also induce changes in calmodulin's conformational landscape along coordinates not probed by vibrational spectroscopy, remaining invisible without additional perturbation of binding site structure. Comparison of FTIR line shapes with 2D IR diagonal slices provides a clear example of how nonlinear spectroscopy produces well-resolved line shapes, refining otherwise featureless spectral envelopes into more informative vibrational spectra of proteins.

A Novel Peptide Restricts Ethanol Modulation of the BK Channel In Vitro and In Vivo.

Alcohol is a widely used and abused substance. A major unresolved issue in the alcohol research field is determining which of the many alcohol target proteins identified to date is responsible for shaping each specific alcohol-related behavior. The large-conductance, calcium- and voltage-activated potassium channel (BK channel) is a conserved target of ethanol. Genetic manipulation of the highly conserved BKα channel influences alcohol-related behaviors across phylogenetically diverse species that include worm, fly, mouse, and man. A pharmacological tool that prevents alcohol's action at a single target, like the BK channel, would complement genetic approaches in the quest to define the behavioral consequences of alcohol at each target. To identify agents that specifically modulate the action of ethanol at the BK channel, we executed a high-throughput phagemid-display screen in combination with a Caenorhabditis elegans behavioral genetics assay. This screen selected a novel nonapeptide, LS10, which moderated acute ethanol intoxication in a BK channel-humanized C. elegans strain without altering basal behavior. LS10's action in vivo was dependent upon BK channel functional activity. Single-channel electrophysiological recordings in vitro showed that preincubation with a submicromolar concentration of LS10 restricted ethanol-induced changes in human BKα channel gating. In contrast, no substantial changes in basal human BKα channel function were observed after LS10 application. The results obtained with the LS10 peptide provide proof-of-concept evidence that a combined phagemid-display/behavioral genetics screening approach can provide novel tools for understanding the action of alcohol at the BK channel and how this, in turn, exerts influence over central nervous system function.

Synthetic Analogues of the Snail Toxin 6-Bromo-2-mercaptotryptamine Dimer (BrMT) Reveal That Lipid Bilayer Perturbation Does Not Underlie Its Modulation of Voltage-Gated Potassium Channels.

Drugs do not act solely by canonical ligand-receptor binding interactions. Amphiphilic drugs partition into membranes, thereby perturbing bulk lipid bilayer properties and possibly altering the function of membrane proteins. Distinguishing membrane perturbation from more direct protein-ligand interactions is an ongoing challenge in chemical biology. Herein, we present one strategy for doing so, using dimeric 6-bromo-2-mercaptotryptamine (BrMT) and synthetic analogues. BrMT is a chemically unstable marine snail toxin that has unique effects on voltage-gated K+ channel proteins, making it an attractive medicinal chemistry lead. BrMT is amphiphilic and perturbs lipid bilayers, raising the question of whether its action against K+ channels is merely a manifestation of membrane perturbation. To determine whether medicinal chemistry approaches to improve BrMT might be viable, we synthesized BrMT and 11 analogues and determined their activities in parallel assays measuring K+ channel activity and lipid bilayer properties. Structure-activity relationships were determined for modulation of the Kv1.4 channel, bilayer partitioning, and bilayer perturbation. Neither membrane partitioning nor bilayer perturbation correlates with K+ channel modulation. We conclude that BrMT's membrane interactions are not critical for its inhibition of Kv1.4 activation. Further, we found that alkyl or ether linkages can replace the chemically labile disulfide bond in the BrMT pharmacophore, and we identified additional regions of the scaffold that are amenable to chemical modification. Our work demonstrates a strategy for determining if drugs act by specific interactions or bilayer-dependent mechanisms, and chemically stable modulators of Kv1 channels are reported.

Coordination to lanthanide ions distorts binding site conformation in calmodulin.

The Ca2+-sensing protein calmodulin (CaM) is a popular model of biological ion binding since it is both experimentally tractable and essential to survival in all eukaryotic cells. CaM modulates hundreds of target proteins and is sensitive to complex patterns of Ca2+ exposure, indicating that it functions as a sophisticated dynamic transducer rather than a simple on/off switch. Many details of this transduction function are not well understood. Fourier transform infrared (FTIR) spectroscopy, ultrafast 2D infrared (2D IR) spectroscopy, and electronic structure calculations were used to probe interactions between bound metal ions (Ca2+ and several trivalent lanthanide ions) and the carboxylate groups in CaM's EF-hand ion-coordinating sites. Since Tb3+ is commonly used as a luminescent Ca2+ analog in studies of protein-ion binding, it is important to characterize distinctions between the coordination of Ca2+ and the lanthanides in CaM. Although functional assays indicate that Tb3+ fully activates many Ca2+-dependent proteins, our FTIR spectra indicate that Tb3+, La3+, and Lu3+ disrupt the bidentate coordination geometry characteristic of the CaM binding sites' strongly conserved position 12 glutamate residue. The 2D IR spectra indicate that, relative to the Ca2+-bound form, lanthanide-bound CaM exhibits greater conformational flexibility and larger structural fluctuations within its binding sites. Time-dependent 2D IR lineshapes indicate that binding sites in Ca2+-CaM occupy well-defined configurations, whereas binding sites in lanthanide-bound-CaM are more disordered. Overall, the results show that binding to lanthanide ions significantly alters the conformation and dynamics of CaM's binding sites.

International Union of Basic and Clinical Pharmacology. C. Nomenclature and Properties of Calcium-Activated and Sodium-Activated Potassium Channels.

A subset of potassium channels is regulated primarily by changes in the cytoplasmic concentration of ions, including calcium, sodium, chloride, and protons. The eight members of this subfamily were originally all designated as calcium-activated channels. More recent studies have clarified the gating mechanisms for these channels and have documented that not all members are sensitive to calcium. This article describes the molecular relationships between these channels and provides an introduction to their functional properties. It also introduces a new nomenclature that differentiates between calcium- and sodium-activated potassium channels.

Structural identifiability of equilibrium ligand-binding parameters.

Understanding the interactions of proteins with their ligands requires knowledge of molecular properties, such as binding site affinities and the effects that binding at one site exerts on binding at other sites (cooperativity). These properties cannot be measured directly and are usually estimated by fitting binding data with models that contain these quantities as parameters. In this study, we present a general method for answering the critical question of whether these parameters are identifiable (i.e., whether their estimates are accurate and unique). In cases in which parameter estimates are not unique, our analysis provides insight into the fundamental causes of nonidentifiability. This approach can thus serve as a guide for the proper design and analysis of protein-ligand binding experiments. We show that the equilibrium total binding relation can be reduced to a conserved mathematical form for all models composed solely of bimolecular association reactions and to a related, conserved form for all models composed of arbitrary combinations of binding and conformational equilibria. This canonical mathematical structure implies a universal parameterization of the binding relation that is consistent with virtually any physically reasonable binding model, for proteins with any number of binding sites. Matrix algebraic methods are used to prove that these universal parameter sets are structurally identifiable (SI; i.e., identifiable under conditions of noiseless data). A general approach for assessing and understanding the factors governing practical identifiability (i.e., the identifiability under conditions of real, noisy data) of these SI parameter sets is presented in the companion paper by Middendorf and Aldrich (2017. J. Gen. Physiol. https://doi.org/10.1085/jgp.201611703).

The structure of binding curves and practical identifiability of equilibrium ligand-binding parameters.

A critical but often overlooked question in the study of ligands binding to proteins is whether the parameters obtained from analyzing binding data are practically identifiable (PI), i.e., whether the estimates obtained from fitting models to noisy data are accurate and unique. Here we report a general approach to assess and understand binding parameter identifiability, which provides a toolkit to assist experimentalists in the design of binding studies and in the analysis of binding data. The partial fraction (PF) expansion technique is used to decompose binding curves for proteins with n ligand-binding sites exactly and uniquely into n components, each of which has the form of a one-site binding curve. The association constants of the PF component curves, being the roots of an n-th order polynomial, may be real or complex. We demonstrate a fundamental connection between binding parameter identifiability and the nature of these one-site association constants: all binding parameters are identifiable if the constants are all real and distinct; otherwise, at least some of the parameters are not identifiable. The theory is used to construct identifiability maps from which the practical identifiability of binding parameters for any two-, three-, or four-site binding curve can be assessed. Instructions for extending the method to generate identifiability maps for proteins with more than four binding sites are also given. Further analysis of the identifiability maps leads to the simple rule that the maximum number of structurally identifiable binding parameters (shown in the previous paper to be equal to n) will also be PI only if the binding curve line shape contains n resolved components.

Closed state-coupled C-type inactivation in BK channels.

Ion channels regulate ion flow by opening and closing their pore gates. K(+) channels commonly possess two pore gates, one at the intracellular end for fast channel activation/deactivation and the other at the selectivity filter for slow C-type inactivation/recovery. The large-conductance calcium-activated potassium (BK) channel lacks a classic intracellular bundle-crossing activation gate and normally show no C-type inactivation. We hypothesized that the BK channel's activation gate may spatially overlap or coexist with the C-type inactivation gate at or near the selectivity filter. We induced C-type inactivation in BK channels and studied the relationship between activation/deactivation and C-type inactivation/recovery. We observed prominent slow C-type inactivation/recovery in BK channels by an extreme low concentration of extracellular K(+) together with a Y294E/K/Q/S or Y279F mutation whose equivalent in Shaker channels (T449E/K/D/Q/S or W434F) caused a greatly accelerated rate of C-type inactivation or constitutive C-inactivation. C-type inactivation in most K(+) channels occurs upon sustained membrane depolarization or channel opening and then recovers during hyperpolarized membrane potentials or channel closure. However, we found that the BK channel C-type inactivation occurred during hyperpolarized membrane potentials or with decreased intracellular calcium ([Ca(2+)]i) and recovered with depolarized membrane potentials or elevated [Ca(2+)]i Constitutively open mutation prevented BK channels from C-type inactivation. We concluded that BK channel C-type inactivation is closed state-dependent and that its extents and rates inversely correlate with channel-open probability. Because C-type inactivation can involve multiple conformational changes at the selectivity filter, we propose that the BK channel's normal closing may represent an early conformational stage of C-type inactivation.

Conserved properties of individual Ca2+-binding sites in calmodulin.

Calmodulin (CaM) is a Ca(2+)-sensing protein that is highly conserved and ubiquitous in eukaryotes. In humans it is a locus of life-threatening cardiomyopathies. The primary function of CaM is to transduce Ca(2+) concentration into cellular signals by binding to a wide range of target proteins in a Ca(2+)-dependent manner. We do not fully understand how CaM performs its role as a high-fidelity signal transducer for more than 300 target proteins, but diversity among its four Ca(2+)-binding sites, called EF-hands, may contribute to CaM's functional versatility. We therefore looked at the conservation of CaM sequences over deep evolutionary time, focusing primarily on the four EF-hand motifs. Expanding on previous work, we found that CaM evolves slowly but that its evolutionary rate is substantially faster in fungi. We also found that the four EF-hands have distinguishing biophysical and structural properties that span eukaryotes. These results suggest that all eukaryotes require CaM to decode Ca(2+) signals using four specialized EF-hands, each with specific, conserved traits. In addition, we provide an extensive map of sites associated with target proteins and with human disease and correlate these with evolutionary sequence diversity. Our comprehensive evolutionary analysis provides a basis for understanding the sequence space associated with CaM function and should help guide future work on the relationship between structure, function, and disease.

Na+ channel function, regulation, structure, trafficking and sequestration.

This paper is the second of a series of three reviews published in this issue resulting from the University of California Davis Cardiovascular Symposium 2014: Systems approach to understanding cardiac excitation-contraction coupling and arrhythmias: Na(+) channel and Na(+) transport. The goal of the symposium was to bring together experts in the field to discuss points of consensus and controversy on the topic of sodium in the heart. The present review focuses on Na(+) channel function and regulation, Na(+) channel structure and function, and Na(+) channel trafficking, sequestration and complexing.

Analyzing single-molecule time series via nonparametric Bayesian inference.

The ability to measure the properties of proteins at the single-molecule level offers an unparalleled glimpse into biological systems at the molecular scale. The interpretation of single-molecule time series has often been rooted in statistical mechanics and the theory of Markov processes. While existing analysis methods have been useful, they are not without significant limitations including problems of model selection and parameter nonidentifiability. To address these challenges, we introduce the use of nonparametric Bayesian inference for the analysis of single-molecule time series. These methods provide a flexible way to extract structure from data instead of assuming models beforehand. We demonstrate these methods with applications to several diverse settings in single-molecule biophysics. This approach provides a well-constrained and rigorously grounded method for determining the number of biophysical states underlying single-molecule data.

Determination of parameter identifiability in nonlinear biophysical models: A Bayesian approach.

A major goal of biophysics is to understand the physical mechanisms of biological molecules and systems. Mechanistic models are evaluated based on their ability to explain carefully controlled experiments. By fitting models to data, biophysical parameters that cannot be measured directly can be estimated from experimentation. However, it might be the case that many different combinations of model parameters can explain the observations equally well. In these cases, the model parameters are not identifiable: the experimentation has not provided sufficient constraining power to enable unique estimation of their true values. We demonstrate that this pitfall is present even in simple biophysical models. We investigate the underlying causes of parameter non-identifiability and discuss straightforward methods for determining when parameters of simple models can be inferred accurately. However, for models of even modest complexity, more general tools are required to diagnose parameter non-identifiability. We present a method based in Bayesian inference that can be used to establish the reliability of parameter estimates, as well as yield accurate quantification of parameter confidence.